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Abstract 

A methodology is presented for investigating the dynamics of heterogeneous media using the nonlocal 
continuum model given by the peridynamic formulation. The approach presented here provides the ability 
to model the macroscopic dynamics while at the same time resolving the dynamics at the length scales 
of the microstructure. Central to the methodology is a novel two-scale evolution equation. The rescaled 
solution of this equation is shown to provide a strong approximation to the actual deformation inside the 
peridynamic material. The two scale evolution can be split into a microscopic component tracking the 
dynamics at the length scale of the heterogeneities and a macroscopic component tracking the volume 
averaged (homogenized) dynamics. The interplay between the microscopic and macroscopic dynamics 
is given by a coupled system of evolution equations. The equations show that the forces generated by 
the homogenized deformation inside the medium are related to the homogenized deformation through a 
history dependent constitutive relation. 

1 Introduction 

The peridynamic formulation introduced in Silling [2A\ is a non-local continuum theory for deformable 
bodies. Material particles interact through a pairwise force field that acts within a prescribed horizon. 
Interactions depend only on the difference in the displacement of material points and spatial derivatives 
in the displacement are avoided. This feature makes it an attractive model for the autonomous evolution 
of discontinuities in the displacement for problems that involve cracks, interfaces, and other defects, see 
[2 m [m [25l [26j |2I] ■ Recent investigations aimed toward developing the numerical implementation, and 
application areas of the peridynamic model include [7] , [30] , [21] , [32] ; [33] ■ More mathematically related 
investigations address issues related to the function space setting of peridynamics [12], [H] and the link 
between the linearized peridynamic formulation and the operators appearing in the Navier system of linear 
elasticity in the limit of vanishing non- locality [12], [29]. In this context the convergence of the solutions 
of the peridynamic equations to the solutions of the Navier system is demonstrated in 8J. In other related 
work the development of a non-local vector calculus with applications to non-local boundary value problems 
has been carried out in 16J. Recent work on the multi-scale applications of peridynamics have shown how 
the peridynamic equations formulated at mezo-scales can be recovered by a suitable upscaling of atomistic 
formulations, see [23] . 

In this paper new tools are developed for the analysis of heterogeneous peridynamic media involving two 
distinct length scales over which different types of peridynamic forces interact. The setting treated here 
involves a long range peridynamic force law perturbed in space by an oscillating short range peridynamic 
force. The oscillating short range force represents the presence of heterogeneities. It is also assumed that 
there is a sharp density variation associated with the heterogeneities. In this treatment we carry out the 
analysis in the small deformation setting. For this case the reference and deformed configurations are taken 
to be the same and both long and short range forces are given by linearizations of the peridynamic bond 
stretch model introduced in ,24j . 
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The relative length scale over which the short range forces interact is denoted by e and points inside the 
domain containing the heterogeneous material are specified by x. Here we will suppose the heterogeneities 
are periodically dispersed on the length scale £ — for some choice of rt = 1, 2, . . . The deformation inside 
the medium is both a function of space and time t and is written The multi-scale analysis of the 

peridynamic formulation proceeds using the concept of two-scale convergence, introduced and developed by 
Nguetseng [20] and subsequently in Allaire [T] , see also E [5] . The two-scale convergence originally introduced 
in the context of partial differential equations turns out to provide a natural setting for identifying both the 
coarse scale and fine scale dynamics inside peridynamic composites. The theory and application of the two- 
scale convergence is taken up in section three of this paper where a novel two-scale peridynamic equation is 
derived. The two-scale formulation is described by introducing a rescaled or microscopic variable y = x/e. 
The solution of the two-scale dynamics is a deformation u{x,y,t) that depends on both variables x and y. 

The rescaled solution u{x,x/e,t) is shown to provide a strong approximation to the actual deformation 



u'^{x,t) inside the peridynamic material. This is is shown in section 3.3 where an evolution law for the error 
e'^{x,t) = u'^{x,t) — u(x,x/£,t) is developed. It is shown that e^{x,t) vanishes in the norm, with respect to 
the spatial variables, when the length scale of the oscillation tends to zero for all p in the interval 1 < p < oo. 
The advantage of using the two-scale dynamics as a computational model is that it has the potential to lower 
computational costs associated with the explicit peridynamic modeling of millions of heterogeneities. This 
issue is discussed in section [3?3l 

It is important for the modeling to recover the dynamics that can be measured by strain gages or other 
macroscopic measuring devices. Typical measured quantities involve averages of the deformation u^{x,t) 
taken over a prescribed region V with volume denoted by To this end we denote the unit period cell for 
the heterogeneities by Y and project out the fiuctuations by averaging over y and write 

u"{x,t)^ / u{x,y,t)dy. (1.1) 



Y 



In section four it is shown that 



lim / u^x,t)dx = -r^ I u"{x,t)dx. (1.2) 



\V\ Jv ' ' ' \V\ 



V 



In this way we see that the average deformation is characterized by u^{x,t) when the scale e of the mi- 
crostructure is small. We split the deformation into microscopic and macroscopic parts and write u{x, y, t) = 
u^{x, t) + r{x, y, t). The interplay between the microscopic and macroscopic dynamics is given by a coupled 
system of evolution equations for and r. The equations show that forces generated by the homogenized 
deformation inside the medium are related to the homogenized deformation through a history dependent 
constitutive relation. The explicit form of the constitutive relation is presented in section four where we 
present a homogenized evolution equation for the coarse scale dynamics written exclusively in terms of , 
see dilTl). 
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Figure 2: Deformation of a bond within the peridynamic horizon. 



1.1 Peridynamic Formulation of Continuum Mechanics in Heterogeneous Me- 
dia 

We consider elastic deformations inside a body described by the bounded domain fi. In the peridynamic 
theory, the time evolution of the displacement vector field u, in a homogeneous body of constant density p 
is given by the partial intcgro-differential equation 

pdfu{x,t)= / f{u{x,t) — u{x,t),x — x,x)dx + b{u,x,t), {x,t) ^ Q x {0,T) (1-3) 

JH-,{x)nn 

where H^{x) is a neighborhood of x of diameter 27, 6 is a prescribed loading force density field, and is a 
bounded set in . Here / denotes the pairwise force field whose value is the force vector (per unit volume 
squared) that the particle at x exerts on the particle at x. For a homogeneous medium / is of the form 
f{u{x,t) — u{x,t),x — x), i.e., it depends only on the relative position of the two particles. We will often 



refer to / as a bond force. Only points x inside H^ix) interact with x. Equation (1.3) is supplemented with 



initial conditions for u{x^ 0) and dtu{x, 0). For the purposes of discussion it will be convenient to set 

£,^x~x, 

which represents the relative position of these two particles in the reference configuration, and 

r] — u{x, t) ~ u{x, t), 

which represents their relative displacement (see Figure [2]). In this treatment, all elastic deformations are 
assumed small and the reference and deformed configurations are taken to be the same. 

We now introduce the heterogeneous peridynamic material. One can think of it as a material with long 
range peridynamic forces acting over a neighborhood of diameter 27 perturbed by an oscillating density 
fluctuation and oscillatory short range bond force acting over a much smaller neighborhood of diameter 2e5. 
Both the long and short range pairwise elastic forces will be given by the linearized version of the bond-stretch 
model proposed in '27] . The long range force is given by 

I 0, otherwise. 

Here ^ (8) ^ is a rank one matrix with elements (^ ® = and 7 is the prescribed peridynamic horizon 
and A is a positive constant. 

In this paper we assume that oscillations in the density and short range bond force are periodic. Here 
the oscillations are characterized by rescalings of a unit periodic peridynamic bond force and density. To 
describe these we introduce the unit period cube F C M'^ for the microstructure. The local coordinates 
inside Y are denoted by y with the origin at the center of the unit cube. The unit cube is composed of two 
or more peridynamic materials with different densities. To fix ideas one can consider reinforced composites 
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made up of an inclusion phase such as a particle or fiber and a second host phase that surrounds the particle 
or fiber. A fiber reinforced material is portrayed in Figure [4] The presence of material heterogeneity is 
reflected by the appearance of peridynamic forces acting within the length scale of the period. Let Xf denote 
the indicator function of the set occupied by the inclusion material and Xm denote the the indicator function 
of the set occupied by the host or matrix material. Here Xf is given by 



1, y is in the inclusion phase, 
0, otherwise, 



Xf{y) = 

and Xm is given by 

Xm(y) = 1 - Xf{y)- 

We extend the functions Xf and Xm to by periodicity. For future reference, we denote by 6{ and the 
volume fractions of the included material and the matrix material, respectively. Here 9{ = Jy Xf{y)dy and 
= 1 — 6i. The density of the matrix material inside the unit period cell is given by the unperturbed 
density = p and that of the inclusion is given by < = p + Ap where Ap can be a positive or negative 
constant. The density characterizing the heterogeneous medium is 

pi.y) = Xf{y)pf + Xm{y)Pm- (i-4) 

The short-range pairwise force is characterized by a bond strength as associated with a horizon 5 > 0. The 
peridynamic horizon 5 is chosen to be smaller than the spacing separating the inclusions. In addition the 
inclusions are assumed to be sufiiciently smooth so that the points y and y are separated by at most one 
interface when Ij/ — < 5. For any two points y and y in M'^ as is given by 



asiy,y) 



Cf, if y and y are in the same inclusion and |?/ — y| < 6 

Cm, if y and y are in the matrix phase and |y — ^| < S 

Ci, if y and y are separated by an interface and \y — y\ < S 

0, if\y-y\>S. 



The material parameters Cf and Cm are intrinsic to each phase and can be determined through experiments. 
Bonds connecting particles in the different materials are characterized by Ci, which can be chosen such that 
Cf > Ci > Cm > 0, see [17]. Mathematically we express the bond strength as 

as{y,y) = xs{y - y)a{y,y), (1-5) 

where xs{z) — 1 for \z\ < S and xsi^) = for \z\ > 5 and a is given by 

a{y, y) = Cf xi{y)xi{y) + Cm Xm{y)xm{y) + G {xi{y)xm{y) + Xm{y)xf{y)) ■ (i-6) 

The short-range peridynamic force deflned on fl is given by 

fshortiv, ^, ^ —a^s j -^V- (1-7) 



For future reference we see from (1.5) and (1.6) that a^si^^ f ) is given by 

aes (^|, = Xes{x - x) (Cf xl{x)xl{i) + Cm xL{x)xL{i) + Ci {xnx)xL{x) + xl{x)xl{i))) , (1-8) 

where Xf (^) •= Xf(f ) and Xm(^) ■= Xm(f )■ The oscillating density for the heterogeneous medium is 
given by pe{x) = p(f ). 

The elastic displacement inside the heterogeneous body f2 is denoted by u^(x,t) and the peridynamic 
equation of motion for the heterogeneous medium is given by 

Pe{x)dfu''{x,t) = / fiong{u%X,t) - U''{x,t),^)dx 



fLoAiu^x,t) -u%x,t)),^,x)dx ^^-^^ 



+ 6^ {x, t) , for X in fi. 
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The peridynamic equation is supplemented with initial conditions 

u%x,0) = ul{x) (1.10) 
dtu'{x,0) = v'oix). (1.11) 

Here the body force b^{x, t) and initial conditions Uq{x), Vq{x) can depend upon e. When these functions are 
bounded in LP(ri)^ for p > 1 it follows from the theory of semigroups that there is a cl assic solution u'^{x, t) 
belonging to C^([0,r]; L'^iVt)'^). This is discussed in the following section, see Remark 



2.2 



In what follows we will develop strong approximations for solutions when the prescribed body forces 
and initial conditions are continuous at the coarse length scale but possess discontinuous oscillations over fine 
length scales. For this choice we look for a solution u'^(x,t) continuous in time but possibly discontinuous 
in the spacial variables and belonging to the Lebesgue space LP{n)^ for 1 < p < oo. In this paper we show 
that we can find solutions u'^{x,t) and strong approximations of the form u{x,x/e,t) that both belong to 
C'^{[0,T]; LP{n)^), for a wide class of initial conditions and body forces. In order to describe this class of 
initial conditions and body forces we consider the space Lp^j.{Y; C{il)^) of functions il^{x, y) measurable with 
respect to y, L^'-integrable on Y and F-periodic in y, with values in the Banach space C{VlY of continuous 
vector fields on Q.. Every element 'ip{x,y) of this space is a Caratheodory function and hence i^ix, |) is 
measurable on Q. and belongs to L'p{{1). This kind of function space is well known in the context of two-scale 
convergence see, [1], and [21]. In what follows we will suppose b{x,y,t) belongs to C([0,T]; LP^j.(Y;C{n)^) 
and both uo{x,y) and VQ{x,y) belong to Lp^j.{Y;C{U)^). For this choice the initial conditions and body 
forces are given by u^{x, 0) — uo{x, '^), dtu'^{x, 0) = vo{x,^), and b'^{x, t) = b{ x, f , t). The construction of a 
strong approximation for this class of data is given in Theorem |3.12| of section |3.3[ 

It is important at this stage to point out that it is precisely the scaling of the bond force together 



with the scaling eS of the horizon that ultimately delivers the macroscopic equations for given by (4.171 
In this context we expect other types of macroscopic equations to arise for different scalings of the bond force 
strength. Recent work for homogeneous media show that the classical equations of linear elasticity arise for 
bond force scaling on the order of and horizons with scaling e, see [H], and [S]. 

When the initial conditions and body force are continuous functions and the density and bond forces 
characterized by Q;(|, |) are also continuous then the solution u'^ is continuous in space and belongs to 
C^([0,T]; C{^)^); this is discussed in the next section. 

In forthcoming work we will focus on the development of strong approximations for initial conditions 
that are discontinuous with respect to coarse length scales. This will be carried out for heterogeneous 
peridynamic media characterized by oscillatory but continuous densities and bond forces. More generally 
one could contemplate strong approximations for more general combinations of bond forces and initial data. 
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(a) (b) 
Figure 4: (a) Composite cube Y. (b) Cross-section of Y along the fiber direction. 

2 Peridynamic Formulation for Heterogeneous Media: A Well 
Posed Problem 

In this section, we make use of the semigroup theory of operators to show the existence and uniqueness of 
solutions to ([L9|-( |1.11| ). For v e LP{n)^, with 1 < p < oo, let 

A^M^) = X^^^^l^^^p^vix)dx, (2.1) 
JH^{x)nn \x — x\ 

AUv{x) = p:\x)[ X^^^^^l^^^f^dxvixl (2.2) 
JH^{x)nn \x — x\ 

H^six)nn £^ \ e e 

H^s{x)nn \s e 

A| = A%,~Ay, (2.6) 
A" = Al+A'g. (2.7) 

Then by making the identifications u^{t) = u^{-, t) and b^{t) — 6(-, t), we can write ( |1.9[ )-( [m| ) as an operator 
equation in L^^n)^ 

u^{t) = A''u^{t)+ p-^b^{t), t€[0,r] 

u'{Q) = ul, (2.8) 

or equivalently, as an inhomogeneous Abstract Cauchy Problem in 
LP{n)^ X LPiflf 

il^t) = A^U^t) + B^t), te[0,T] 

where 



AF / N -1/ \ ; 1 f X x\ (x - x) ® (x — x) . ^. ,^ ^. 

A%Mx) = p,\x) I -a( -,-)'- ^^^^ '-v{x)dx, (2.3) 

AF / \ -1 / \ ; 1 f X x\ (x - x) ^ (x — x) . . ^, 

^l,2^'(a^) = Pe i^) I 72'^{-^-]- u:^-^3 -dxv{x). (2.4) 



Also we set 



(2.9) 



- ( (') ) ' ^0 n 1) ' = ( ,-.2nt) ) ' - ( i 

Here / denotes the identity map in L^(il)'^. 

Proposition 2.1. Letl<p<oo and assume that b G C([0, T]; ^^^^(F; C{Uf)) and Uo E i^er(^; C"!^)'^) x 
L^,^(y; (7(0)3). y^jg^ 
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(a) The operators A'^ and M are linear and hounded on LP(f2)^ and LP{n)^ x LP{n)^, respectively. More- 
over, the bounds are uniform in e. 



(h) Equation (2.9) has a unique classical solution in C^([0,T]; LP{VL)'^ x U'{Q.)'^) which is given by 

U'{t)^e'^'U^+ f e^'^''^^' B%t) dT, ie[0,T], 
Jo 



(2.10) 



whe 



(2.11) 



Moreover, equation (2.8) has a unique classical solution G C^([0,T]; L^(r2)'^) which is given by 



u' (t) = cosh (^t\^^ Uq + VA" ^ sinh (tV^^ 



with the notation 



A^ / sinh {{t-T)VAnb%T)dT 



(2.12a) 

(2.12b) 
(2.12c) 



Remark 2.2. The hypothesis of Proposition \2.1\ can be relaxed by assuming that the sequences of initial 
conditions (mq),(wo), are bounded in U'(Q.Y o.'^d (&^(-,i)) is uniformly hounded in L'''{Vl)'^ for t G [0,r]. 
This is proved following the same steps given in the proof of Proposition \2.1\ presented below. 



n=0 ^ ' 
-1 / \ °° f2n+l 

/A? sinh ( y (A^)" 



(c) The sequences (u'^)e>Oi {'^'^)e>o, O'f^d {u'^)e>o o,'"'^ hounded in _L°°([0,r]; L^iJlY). 



Proof. Part (a). It is clear that the operators ^, ^S2' ^1 1' 2 ^-re Unear. So we begin the proof 

by showing that A'^g and A'y ^ ai'c uniformly bounded sequences of operators on LP[^^ for 1 < p < 00. We 
introduce the indicator function XiA^) taking the value one for x inside f2 and zero for x outside il and let 
V denote a generic vector field belonging to U'i^Y ■ Then by the change of variables x — x + ez vcl (2.3) we 
obtain 



A%^v{x) = / Xfi(a; + £2)a f-, - + ^— i;(a; + £z) dz. 

Applying Minkowski's inequality gives 



(2.13) 



i/p 

v{x + ez)\Pdx] dz. (2.14) 



Let a — max p ^{y)a{y,y') and we see that 
v.v'eY 



P!s,i«(a^)IUp(o)3 < a 



1 



Hs{0) Pl 

< Ms\\v\\LP(nr, 



i/p 

Xn{x + ez)\v{x + ez)\P dx \ dz 



where M5 is independent of e and given by 

M5 = a f 



Hs(0) 



dz \ — a 



27r<52 



(2.15) 



(2.16) 
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which shows that the operators are Ag ^ is uniformly bounded with respect to e. Similarly, Ag 2 can be 
written as 

Al^^vix) = / Xn{x + ez)p-\-)a(-,-+z]^r^dz v{x). (2.17) 



HsiO) 



Thus 

\Ayv{x)\ < Ms\v{x)\, 

from which the boundedness of Ag^ immediately follows. Combining these results shows that Ag, which is 
given by Ag — Ag 2, is a sequence of uniformly bounded operators on LP{U)^. 

Next we show that the linear operators A^j^ — A'j^ ^ — A'j^ ^ are a sequence of uniformly bounded operators 
on LPi^Y ■ Changing variables x = x + ^ and applying Minkowski's inequality gives 

Jh.,{o) \Jn e 141 / 

< ML\\v\\LP(nr (2-18) 



where is given by 



Ml = ma^{p-'iy)} x A^, (2.19) 



and it follows that the operator A"]^ is bounded in LP{n)^. The boundedness of A"]^ 2> which is given by 
(2.2), follows immediately from its definition. Therefore is imiformly bomided on LP{il)^ with respect 
to £. 

Since A'^ — A'^ + Ag , we conclude that 

\\A'v\\LP(nr < M\\v\\L.(nr, (2-20) 

for a positive constant AI which is independent of e. The operator is clearly linear, thus it remains to show 
that this operator is uniformly bounded on LP{Q)'^ x LP{fl)^. To see this, we let {v,w) S LP{n)^ x LP{fl)^. 
The norm in this Banach space is given by 



\\iv,w)\\LP{n)3xLp{n)3 = ||w||Lp(n)3 + \\w\\Lp{n)3- 



We note that 



Thus we obtain 



A" 



V \ _ / I \ / V \ _ / w 
w I ~ [ I [ w I ^ [ A'v 



||A''(w, w)|jLP(o)3xLp(o)3 = \\w\\LP{n)3 + \\A%\\Lp(^n)3 



< 



kllLp(a)3 + llA-|llkl|Lp(o)3. (2.21) 



From ( |2.21[ ) it follows that 

w)||LP(f2)3xLP(f2)3 < w)||LP(n)3xLP(n)3, (2.22) 

for some positive constant M completing the argument. 

Part (b). We have seen from Part (a) that A*^ is a bounded linear operator on the Banach space LP{n)'^ x 
LP{n)^. Also, since b" is in C([0,T]; LP{nf), it follows that = (0,6-^) is in C{[0,T]; LP{nf x LP{n)^). 
From these facts, it follows from the theory of semigroups, see for example, [T5] . 

1. The operator A^ gene rates a uniformly continuous semigroup {e**^}(>o on LP(SI)'^ x L^(r2)'^, where 
gtA^ is given by (2.11 1. 
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2. The inhomogeneous Abstra ct Ca uchy Problem ( 2.9 1 has a unique classical solution £ C^{[0,T]; LP{il)^x 
LP{n)^) which is given by ( |2.10[ ). 

It immediately follows from ^ that the second order inhomo geneo us Abstract Cauchy Problem (2.8) has a 
unique classical solution u'^ E C^([0,T]; LP{n)^) and formula (2.12) follows immediately from (2.11). 
Part (c). We recall that 

X 

Uq{x) := uo(x, -) 
e 

X 

v^{x) := Vo{x,-) 
e 

where uo{x, y), vo{x, y) are in LP^^{Y; C{fl)^). We surround by a cube of integer side length L and extend 
uo{x, y) to the cube by setting UQ{x,y) — for x outside Q and for every yinY. We note that the extended 



is e 



periodic in the second variable and shift the cube so that it is commensurate with the 



periods. The period cells of side length e are denoted by eYi and the cube is given by their union UisYi 
where the index i ranges from 1 to L^n^. Since we have extended UQ{x,y) so that it vanishes when x lies 
outside fl one can write 



Hence 



l|uolkp(n)3 = 



|wollif(n)3 < 



\uo{x, dx 



i/p 



(2.23) 



sup \uo{x' , —)\P dx 



i/p 



V / sup \uo{x', -)\Pdx 



1/p 



(2.24) 



Here the last inequality follows from the change of variables y — j:- Thus Uq is uniformly bounded in LP{n)^. 
Similarly is uniformly bounded which implies that Uq is uniformly bounded in LP{il.)^ x LP{il)^. The same 
considerations show that for t S [0,T], that b'^{t) is uniformly bounded in LP{il)^. Since b^{t) is continuous 
in t, it follows that b'^ is uniformly bounded in C([0, T]; LP(il)^), which implies that is uniformly bounded 
in C([0,r]; LP{n)^ x LP{nf). 
Next we note that 



< 
< 



otiiA^II 



tM 



(2.25) 



where in the last inequality we have used the fact that A*^ is uniformly bounded. Taking the norm in both 



sides of (2.10) and by using (2.25), we obtain 



\U%t)\\ 



(2.26) 



for some positive numbers Mi , M2 , and M. This implies that [/^ is uniformly bounded in L°° ( [0, T] ; LP{il)^x 
LP{n)^). Ther efore the sequences (u^)e>o and (M'^)e>o are bounded in L°°{[0,T]; LP{i})^). Finally, it follows 
from equation (2.8) that the sequence (■u^)e>o is bounded in i°°([0,T]; LP{n)^), completing the proof. □ 



It is easily seen that for continuous initial conditions and body forces that the peridynamic solution is 
also continuous in space provided that the bond forces and densities are continuous. To fix ideas we "smooth 
out" the characteristic functions Xf and Xm by mollification. Indeed given any infinitely differential function 
C with compact support on we fix /3 such that < (3 < S and form C^{x) — /3~^C(f )• The mollified 
characteristic functions are given by Xf{^) — iC^ * Xf){x) and Xmi^) ~ (C^ * Xm){x). The replacement of 
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Xi and Xm by their mollified counter parts in (1.4) and (1.6 1 delivers a short range bond force /JjjQj.j(77, ^, x) 
and density p^{x) that are continuous in x. For this case it is easy to see that -j^, Ai^ i, and A]^ 2 are 

linear operators mapping C{^Y i^to itself. A straight forward application of Holder's inequality shows that 
^Si' ^S2' ^L,l^ and Al^2 are bounded and that the operator norms of A^-^, A^g2 are uniformly bounded 
with respect to e. To fix ideas we choose and uq in (7(0) and for 6 in C^([0, T]; C(f2)) and proceeding as 
before we find that the solution u'^ of the peridynamic initial value problem exists is unique and belongs to 



3 Strong Approximation by Two-Scale Functions 

The aim of this section is to build an approximation of u^{x,t) when the period e of the microstructure is 
small. In what follows we show how to systematically identify a function u{x,y,t) that is oscillatory with 
respect to a new "fast" spatial variable y that when rescaled y = - delivers a strong approximation to 
u'^lx, t), i.e., 

X 

lini \\u''{x,t) - u{x, -,t)\\Lp(i^) = 0. (3.1) 

It is shown that the desired function u{x,y,t) is the "two-scale" limit of the sequence {u'^{x,t)} for e — > 0. 
After periodically extending u{x,y,t) in the y variable we find that it satisfies the two-scale peridynamic 
initial-value problem given in theorem |3.10[ In the subsequent sections we apply this fact to show that 
u{x, -,t) provides a strong approximation to u'^{x,t) when e is sufficiently small. 



3.1 Two-Scale Convergence 

To expedite the presentation we list the following useful function spaces 

/C = {i/; e C^(R^ X Y), i){x,y) is ^-periodic in y}, 
J = {i^ e C~(K3 X y X M+), ^l){x,y,t) is Y-periodic in y}, 
- {u;eC([0,r];LP,,(y;C(n)3)}, 

and introduce the definition of two-scale convergence. Let p and p' be two real numbers such that 1 < p < oo 
and l/p + l/p' = 1. 

Definition 3.1 (Two-scale convergence j20l IT]). A sequence {v^) oj Junctions in L^ijl), is said to two-scale 
converge to a limit v £ {Q x Y) if, as £ — >■ 

v^{x)ip(x, — j dx ^ v{x,y)4'{x,y) dxdy (3-2) 



2 



for all ip £ LP (il; Cper{Y)). We will often use ^ v to denote that {v^) two-scale converges to v. 



If the sequence {v^) is bounded in LP{rt) then (il; Cper^X)) can be replaced by K, in Definition (3.1 ) 
(see [21] )• For time-dependent problems one slightly modifies the above two-scale convergence to allow for 
homogenization with a parameter, see [51 [3] • Here the parameter is denoted by t. 

Definition 3.2. A bounded sequence {v^) of functions in LP{D, x (0,r)), is said to two-scale converge to a 
limit V e LP{n X r X (0, T)) if, as e ^ Q 

{x,t)ilj (x, — ,tj dxdt ^ I v{x,y,t)ip{x,y,t) dxdydt (3-3) 

nx(o,T) V e / JnxYx{o,T) 



for all e J. 



Definition 3.1 is motivated by the following compactness result of Nguetseng, see [50] and Allaire [T]. 
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Theorem 3.3. Let (w^) be a bounded sequence in L^{^1). Then there exists a subsequence and a function 
V e LP{fl X Y) such that the subsequence two-scale converges to v. 

A similar two-scale compactness holds for time dependent problems and is stated in the following theorem. 

Theorem 3.4. Let (u^) be a bounded sequence in L^irt x (0,T))'^. Then there exists a subsequence and a 
function v £ L^iTl xYx {0,T))^ such that the subsequence two-scale converges to v. 

The proof of compactness for the time dependent case is essentially the same as the proof of Theorem 



3.3 A slight variation of Theorem 3.4 can be found in [5] and [5]. For future reference we recall the following 



well known results on two-scale convergence that can be found in |21j . 

Proposition 3.5. Let (w^) be a bounded sequence in LP{fl x (0,r))'^ that two-scale converges to v G LP{fl x 
Y X (0,T))3. Then as e ->■ 

v{x, y, t) dy weakly in L'P{Vt x (0, T)f. 

Proposition 3.6. Ifv^{x) converges to v(x) in LP{fl)^ then its two-scale limit is v. 

Last we state two-scale convergence theorems for test functions. 
Proposition 3.7. Ifij^ix^y) belongs to K or LP,^^{Y]C{Q)'^) then ip(x, |) two-scale converges to 'ip{x,y) and 

lhaUix,^)\\l,^^^^ f \4'{x,y)\Pdxdy. (3.4) 

Moreover given any bounded sequence in LP{Q,)^ two-scale converging to v then 



lirn I {x)ip{x, —) dx = / v{x,y)'ip{x,y) dxdy (3.5) 
^ JnxY 



for every test function belonging to LP^^{Y; C(f2)'^). 

Similarly if ip{x^y,t) belongs to J or Cp then ^{x, - ,t) two-scale converges to ijj(x,y,t) and 



lim llV^a;, -,i)lliP(nx(0T))3 = / \i^{x,y,t)\Pdxdydt. (3.6) 



Moreover given any bounded sequence in LP{Q, x (0,r))'^ two-scale converging to v then 

lim / {x,t)ip{x, — ,t) dxdt = / v{x,y,t)'ip{x,y,t)dxdydt (3.7) 

^^'^JOx(0,T) £ JnxYx{0,T) 

for every test function ip belonging to Cp. 



3.2 The Two-Scale Limit Equation 

In this section, we use two-scale convergence to identify the limit of the solution u^{x,t) of ( |1.9[ )-( LTi] ) for 
initial data Uq — uo{x, |), vq = vo{x, |) and body force b^{x, |,i) with uq and vq in LPg^{Y;C{n)'^) and 
b e Cp. For V e LP{n)^, with § < p < oo, let 

ir t \ f {x-x)(^{x~x) 

Kl,iv{x) = / A v{x)dx, (3.8) 

JH-,{x)nn \x--x\-^ 

H-,(x)nQ, 

H^s(x)r\n £^ \£ s 



I \ I X (* - 2;) (g) (i - 
Kl,2V[x) = / A dxv(x), (3.9) 



j^e , \ I 1 f X x\ {X - X) <S) {X - X) 

Ks^Mx) = / 7^ a (7. 7) 1^_^|3 v{x)dx, (3.10) 

j^e f I 1 f X x\ {X - X) (E) {X - X) 

Ksov{x) = / F n dxv{x). (3.11) 



Bacim Alali and Robert Lipton 



12 



Set Kl — Kl i — 2 and Kg — Kg ^ — Kg ^ and the peridynamic equation (1.9) is written 



p{-^d1u'{x,i) = {KL+KI;)u'{x,t) + b(^x,^,t^ 



(3.12) 



We start by noting that the loading force and initial data are in Cp and LP^^.{Y;C{ft)^) respectively and 
from Proposition |3.7| satisfy the following 



b{x,y,t), 
VQ{x,y). 



(3.13a) 
(3.13b) 
(3.13c) 



We note that from Proposition 2.1 ^c) and Theorem 3.4 it follows that, up to some subsequences, u"^ — ^ u, 
ii'^ ^ u*, and il'^ u** , where u, u* , and u** are in LP{Q xYx [0,T])^. We shall see later that u{x,y,t) 
is uniquely determined by an initial value problem. Therefore u is independent of the subsequence, and the 
whole sequence (w*^) two-scale converges to u. 

We start by extending the function u{x, y, t) in the y variable from Y to as a y-periodic function. 
The next task is to identify the dynamics of the periodically extended u{x,y,t). We multiply both sides of 
(3.12) by a test function ?/;(a;, |, t), where i/'(a;, y, t) is y-periodic in y and is such that ^ G C^iM? x y x M)'^, 
and integrate over Q, x M+ 



OxR+ 



d^u'^{x, t) ■ ij} {x,-,tj p{ — ) dxdt 

({Kl + Kl)u%x, t) + b(^x,^,t^y^P (x, |, t) 
After integrating by parts twice, we obtain 

/ u^{x,t) ■ dftlj (x,-,t) p{ — )dxdt 

JnxR+ \ £ J £ 

+ I w^(x,0) •atV'(a;, ^,o) p(^) dx 

\Kl + Kl)u'{x, t) + 6 (a:, ^ , t) ) • ^ (x, p 



dxdt 



dtu^{x, 0) • %l) {x, -7 p(-) dx 



/axB+ 

Passing to the e — limit we obtain 

u{x, y, t) ■ dtil^ix, y, t)p{y) dxdydt 



dxdt 



vo{x, y) ■ ip{x, y, 0)p{y) dxdy 



OxY 



lim 

e->-0 



QxYxI 



uo{x, y) ■ dtipix, y, Q)p{y) dxdy 

{KL+Kl)u'{x,t)-^(x,-,t) 
b{x, y, t) ■ ip{x, y, t) dxdydt 



dxdt 



(3.14) 



We will use the following lemma to compute the limit on the right hand side of (3.14). 
Lemma 3.8. Let w be in LP^^{Y; C(il)'^) with | < p < 00. and define 

{x — x) <S5 (x — x) 



BLw{x,y) 
Bsw{x,y) 



X- 



\x — x\- 



w{x,y')dy' - ■w{x,y) dx, 



a(y,2/) p {w{x,y) - w{x,y)) dy. 

Hs{y) y-y^ 



Ttien as e — > 0, 
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(a) KLu''{x,t) ^ BLu{x,y,t). _ 
Moreover, the operator p^^Bj^ is linear and bounded on LP^j,(Y; C{fl)^). 

(b) K%u%x,t)^Bsu(x,v,t). 

Moreover, the operator p^^Bs is linear and bounded on LP^^.{Y; C{Q,)^). 

The proof of this lemma is provided at the end of this subsection. 



Remark 3.9. Results similar to Lemma 3.8 can be proven for other function spaces as well. The space 
LperO^'j ^i^)'^) statement of this lemma can, for example, be replaced with the function space LP^^{Y] LP(fl)^) 

or by the function space LP{V,; CperiY)^), where | < p < oo m each of these spaces. 



Apphcation of Lemma (|3.8|) gives 



hm 



{KL+KI,)u'{x,t)-^(^x,'^,t^ 



dxdt 



{Bl + Bs)u{x, y, t) ■ ipix, y, t) dxdydt. 

OxYxR+ 



Thus (3.14) becomes 



axYxi 



u{x, y, t) ■ dfipix, y, t)p{y) dxdydt 



ua{x, y) ■ dtip{x, y, 0)p{y) dxdy 



wo(a;, y) ■ 4'{x, y, 0)p{y) dxdy 



OxY 



{{Bl + Bs)u{x, y, t) + 6(x, y, t)) ■ ip{x, y, t) dxdydt 



(3.15) 



OxYxI 



derivatives with respect to t, for ahnost every t, and the initial conditions supplementing (3.15) are given by 

u{x,y,Qi) = uo{x,y), dtu{x,y,0)=vo{x,y). (3.16) 



Thus by integrating by parts twice, equation (3.15) becomes 



L 



f2xYxK+ 



P[v)d'tu{x, y, t) ■ ip{x, y, t) dxdydt 



{{Bl + Bs)u{x, y, t) + b{x, y, t)) ■ i){x, y, t) dxdydt 



(3.17) 



nxYxB+ 



Since this is true for any function -0 G C^(M^ x F x M)'^ for which tl;{x,y,t) is F-periodic in y, we obtain 
that for almost every a;, y, and t 

d^u{x,y,t) = p'\y)Bu{x,y,t) + p-^b{x,y,t), (3.18) 



where B = B^ + Bg. It follows from Lemma 3.8 that p~ B is a bounded linear operator on L'^^^{Y\ C{VLy), 
with I < p < oo. Therefore the initial value problem given by (3.18 ) and (3.16 1, interpreted as a second-order 
inhomogeneous abstract Cauchy problem defined on L^^.^{Y\ C{Vl)^), with body force in Cp, | < p < oo. 
From the theory of semigroups [52 IS] it follows that this problem has a unique solution m(x, y, t) £ Qp, | < 

P < CO. 

The following summarizes the results of this subsection. 



Theorem 3.10. Let {u^) be the sequence of solutions of (L9)-(L11) with initial data Uq = uo{x, 
vo{x, |) and body force &^(a;, |, t) with uq and vq in L'P^^iY; C{^Y) and b G Cp. Then 
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^ u and the periodic extension of u{x,y,t) in the y variable from Y to M"^ also denoted by u belongs to 
Qp, with I < p < oo, and is the unique solution of 

. f ^ ix — x) ® {x — x) f [ , ^ I \ , I / nA,-- 

p{y)dfu{x,y,t) = / A / u{x,y ,t) dy - u{x,y,t) \ dx 

JH-,{x)n^i \x-x\'^ \Jy ) 

a{y, y) — — ^ — — ("(a;, y, t) - u{x, y, t)) dy ^^'^^^ 

+ b{x,y,t), 
supplemented with initial conditions 

u{x,y,0) = uo{x,y), (3.20) 

dtu{x,y,0) ^ VQ{x,y). (3.21) 

We conclude this section by showing that u is twice diffcrentiable with respect to time and proving 
Lemma 13.81 

Lemma 3.11. Let t G [0,T] and define 

g{x,y,t)= f f u**{x,y,l)dldT + tu*ix,y,0)+u{x,y,0). (3.22) 
Jo Jo 

Then g is in LP{il xY x {0,T))^, twice diffcrentiable with respect to t almost everywhere, and satisfies 

(a) For almost every x,y, and t, g{x,y,t) — u{x,y,t), dtg{x,y,t) — u*{x,y,t), 
and dfg{x, y, t) ^u**{x, y, t) . 



(b) For almost every x and y 



g{x, y, 0) = u{x, y, 0) = ua{x, y), 
dt9{x,y,0) = u*{x,y,0) = vn{x,y). 



Proof. Part (a). Let ^i{x,y) be in C^{Q x Y)^ and F-periodic in y, and let (j) be in C^(M+). Then by 
using integration by parts, we see that 



dtu'^{x, t) ■ -01 (x, - ) (f>{t) dxdt = — u^{x, t) ■ tpi (x, - J 0(t) 

^ JnxR+ ^ 

2 

Sending e to and using the fact that, up to a subsequence, dtuF u* , we obtain 

u* (x, y, t) ■ ipi (x, y) (f>{t) dxdydt 

SlxFxR+ 

u{x, y, t) ■ ipi (x, y) (f>{t) dxdydt. 



dxdt. 



lnxYxR+ 

Since this holds for every ipi we conclude that 

u*{x,y,t)(t){t)dt = - I u{x,y,t)4:{t)dt, (3.23) 

for almost every x and y and for every € C^{M.~^). Similarly, by using the fact that, up to a subsequence, 
diU'^ u**, we see that 

u**{x,y,t)<j)(t)dt^ I u{x,y,t)4){t)dt, (3.24) 
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for almost every x and y and for every (j> £ C^(E+). We note that from (3.801 it is easy to see that g 
twice differentiable in t almost everywhere and satisfies 



IS 



dtgix,y,t) = / u**{x,y,T)dT + u*{x,y,0), 
Jo 

dtg{x,y,t) = u**{x,y,t). 



(3.25) 
(3.26) 

We will use these facts together with (3.231 and ( 3.24[ ) to show that dtg — u* almost everywhere and g — u 
almost everywhere. 

For (j)€C^(R+), we integrate by parts using ( |3.26[ ) and ( |3.24[ ) to find that 



R+ 



dtg{x,y,t)4>{t)dt = / u* {x,y,t)(p{t) dt 



R+ 



Thus we obtain 



{dtgix, y, t) ~ u* {x, y, t)) ^{t) dt = 0, 



(3.27) 



for every (j) e C^(M+). Since dtg{x,y,0) = u*{x,y,0), we conclude from (3.27) that dtg{x,y,t) = u*{x,y,t) 
almost everywhere. Finally it easily follows from ( |3.23 1) that 



{g{x, y, t) - u{x, y, t)) (j){t) dt = 0, 



(3.28) 



for every G C^(IR+). Since g{x,y,0) — u{x,y,0), we conclude from (3.28) that g{x,y,t) = u{x,y,t) almost 
everywhere, completing the proof of Part (a). 



that 



■■'--'^"7 jr-- "o — Jr^ v / 

Part (b). Let ip{x,y,t) be in C^{fl xY x M.)^ and F-periodic in y. Then on integrating by parts, we see 



Oxl 



ku%x,t) • V (x, 



dxdt 



OxR+ 



'{x,t) ■ dtip (^x, 



dxdt 



- / u%x,0)-i:{x,^,0^ 



dx. 



Sending e to 0, we obtain 

u* {x,y,t) ■ Tp {x, y, t) dxdydt 



/sixyxR+ 



/ u{x,y ,t) ■ dti' {x,y,t) dxdydt 

JnxYxR+ 



JnxY 



uo{x,y) ■ilj{x,y,0) dxdy. 



(3.29) 



On the other hand, from Part (a), we see that 

u* {x,y,t) ■ ip {x, y, t) dxdydt 



/SlxyxR+ 



/SlxyxR+ 



u{x, y, t) ■ dtip {x, y, t) dxdydt 



JnxY 



u{x,y,0) ■il!{x,y,0) dxdy. 



(3.30) 



From (3.29) and (3.30) we obtain that 

{uo{x, y) - u{x, y, 0)) • V' (2:, y, 0) dxdy = 0, 



jQxY 



for every ip. Therefore 



u{x,y,0) = uo{x,y), 
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almost everywhere. Similarly we can show that 

dtu{x,y,0) = vo{x,y), 

almost everywhere, completing the proof of Part (b). □ 
Proof of Lemma \37^ Part (a). We compute the two-scale limits of K^^iu'^ and i^L,2'«^ to show that as 



K^u^ix, t) Blu{x, y, t). 



(3.31) 



Let V G C;f (E^ X Yf such that il){x,y) is F-periodic in y, and (f) e C^{M+). Then from the definition of 



K^^i, equation (|3.8|), we see that 



axR+ 



dxdt 



A 



{x — x) ® (x — x) 



\x — XV 



u^(a;, t) dx ■ (x, — j (j){t) 



dxdt, 



Since u^(x,t) u{x,y,t), we obtain using Proposition 



3.5 



that, as e — > 0, 



u{x,y,t)dy weakly in LP{fl x {0,T)f 



It follows from (3.331 that, for fixed x, 



lim 



/ . X , — a^) (8) (x — cc) . , „ s , / s , . , 



+ JH.,(x) 



\x — x\ 



(3.32) 



(3.33) 



(3.34) 



Xn{x)X 



{x — x) (x — x) 



u{x, J/', t) dy' ) (j){t) dxdt 



Here xsi is the characteristic function of fi, taking value 1 for x m. and zero outside. Applying Holder's 
inequality for p + = 1 gives 

{x-x)®{x-x) 
Xn{x)\- r:^ u''{x,t)dx 

H.,(x) \X-X['' 



< X 



< A 



Xnix) ^ , dx 
Hs{x) \x~-x\P ^ 

\ i/p' 



i/p' 



i/p 



\u''{x,t)\P dx 



Hs{x) 



1 



— dx 



\L-^{lo,T];Lp(n)^), for almost every t G [0,r] 



(3.35) 



'h,{x) \x-x\p' 

We note that the integral on the right hand side of the last inequality is finite for p' < 3. From Proposition 
I'"^I|l°°([o,t]; LJ'(si)3) is bounded. Thus from (3.341, and (3.351 and by using Lebesgue's dominated 



2.1 



convergence theorem, we conclude that the convergence of the sequence of functions in (3.34) is not only 
point-wise in x convergence but also strong in LP{^1)^ , with | < p < oo. Therefore from Proposition 3.6 and 



(|3.34|) it follows that the limit of (|3.32|) as e — > is given by 

ifL,iu^(a;,i)-v(x,|)0(t) 



lim 



dxdt 



f2xR+ ^ ^' 

Bl,iu{x, y,t) ■ ip (x, y) 4>{t) dxdtdy, 

nxR+xY 



(3.36) 
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where 



BL,iu{x,y,t) 



A 

H-,{x)nn 



{x — x) <^ {x — x) 



u{x, y' , t) dy' I dx 



(3.37) 



depends only on {x,t) and is constant in y. Next we evaluate the two-scale limit of Kl^2U^ ■ We recall from 
(l2?2l) that 



KL^2u''{x,t) 

I H^,(x)nn 

from which immediately follows that as e — > 0, 



{x - x) ® {x - x) 
X- -r- dx u^{x, t), 



X 



— x) (g) (x — x) 



\x — x[ 



dx u{x, y, t) = Bl,2u{x, y, t). 



(3.38) 



(3.39) 



The result (3.31) follows on combining equations (3.361 and (3.391 and writing Bl — Bl^i — -Bl,2. It is 
evident that p~^B]^ is a linear operator on the Banach space L^^^{Y\ C{Vty^). To show boundedness we 



show that p ^Bl.i and p ^Bl^2 are bounded operators on Lp^^{Y; C{fl)^). For w in Lp^j.{Y; C{il)^) we 
write i = a; + and 



(3.40) 



/ P \y) sup 



Xn{x + i)X—— 



w{x + („y')dy' di 



i/p 



dy 



< 



< X 



H-,{0) 



p-\y) 



x_ 



P \ 1/P 



Y JY 



s^v\xn{x + £)w{x + S„y')\dy'd£,\ dy 

Y X£a J y 

P^^iy) snp \w{x,y')\dy] dy' 



(0) 



where the second inequality follows from Minkowski's inequality and it follows that p~^Bl i is bounded. It 
is evident from the definition of i?L.2 that p^^-Bl,2 is a bounded operator on LP^j.(Y; C{ft)^). 

Part (b). Since Kg — Kg ^ — Kg^j we will compute the two-scale limits of Kg -^u^ and Kg ^u'^ , to show 
that as £ — > 0, 



Let 4 i{x,y,t) = ip2{x)il^iiy)4>{t), where '02 G O, 
(2.13), replacing v{x) with u'^{x,t), we have 



K'gu%x,t) ^ Bsu{x,y,t). 

^1 e C^eriYf, and <j> e 



(3.41) 
Then by using 



K 



nxR+ 



liu%x,t) ■ {x, 



dxdt 



, (XX \ z®z , I X \ 

Xn\x -f £z)oL — , — V z] u [X + ez, t) dz ■ [ x, — ,t ] dxdt, 

\e e / \z\-^ \ e J 



(3.42) 



where xn denotes the indicator function of fi. Thus after a change in the order of integration in the right 
hand side of equation (3.42), we see that 



K\ 



nxR+ 



'■yU^{x,t) ■ 4) {x, ^,t^ 



dxdt 



HsiO) 7oxR+ 



Xn(x -I- ez)a (^^,^ + u'^{x + ez,t)-z tJji (J^^-z ■02(a;)0(O 



dxdtdz. 



(3.43) 



Bacim Alali and Robert Lipton 



18 



Now we focus on evaluating the limit as e — > of the inner integral in (3.43). By the change of variables 
r — X + ez we obtain 



Xn 



SlxR+ 



{x + £z)a ^ 



e e 



+ z]v: 



'{x + ez,t)-z ipi --z '02(a;)(/)(i) 



dxdt 



Xn{r)xn{r - ez) a - z, u^{r,t)-z i/'i - zj -z i;2{r - ez)(j){t) 



drdt 



(3.44) 



:=a-(z), 
We will show that for z G i?5(0), 

lima'(z)= / a [y — z,y) u{r,y,t)- z {y — z)- z ip2{r)<p{t) drdydt. 

JnxYxR+ 



(3.45) 



To see this, we approximate xn by smooth functions C„ such that as n — >■ oo, ^„ — >■ hi L^^^(IR^), with 
l/p+ l/p' — 1. Then by adding and subtracting ^„(r — ez) xn{r — ez) in (3.44), we see that 



a^z) 
where, 

n.e ( \ 
Oi (z) 



n.e / \ 

^2 (2) := 



From Proposition [2?T 



Xn(7-) (xo(?- - £z) - Cn(r - £2;)) X 

^ ^) (^i ■ "01 z^ ■ z '>p2{r — ez)(l){t) drdt, 

Xn{r)Cn{r ~ ez) x 

;R+ 

a 1^ z, u^(r, t) -z -i/ii ^ -z -02(?' ~ £2)(/)(t) drdt. 



sup ||u^|lL-([0,T]:LP(n)3) < OO 
e>0 



(3.46) 

(3.47) 

(3.48) 
(3.49) 



So from (3.47) and on application of Holder's inequality, we see for some constants Ci and C2 that 



i/p' 



K'{z)\<Cil^J^Jxnir-ez)-Cnir-ez)\P dr j x \\u'\\L^^^(^M+.,LHn)^) (3-50) 

|a2'^(z)| < C2||u''||l~^(r+;Lp(o)3) (3.51) 

so there is a constant C such that |a^(z)| < C for e > 0. On the other hand, the second factor on the right 
hand side of (3.50) goes to zero uniformly in e as n — cx) and we conclude that for all e > and z G Hs{0), 

lim a"'^(z) = 0. (3.52) 

Now for n fixed we see that as e — >■ 0, Cn(^ ~ ez)'ijj2{r — ez) — > Cn(^)V-'2('') uniformly. Therefore, we see from 
(|3.48|) that 



lim ao''^(z) 
e-i-O ^ 



lim 

£^0 



xn 



(r)C„(r) a 



r r' 
z. - ] 

e ' e. 



{r,t)-z - z^ -z 'ip2{r)(j){t) 



drdt 



Cnir)a {y - z,y) u(r,y,t)-z ipi (y - z)-z ■il)2{r)(j){t) drdydt, 



nxYxR+ 



(3.53) 
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where in the last step the fact that (u'^)e>o two-scale converges to u{r, y, t) was used. By takmg the limit 
n cx) in (3.53), we obtain 



as 



lim lima2''^(^ 



n— >-oo £— ^0 



/i2xyxR+ 



a{y - z, y) u{r, y,t)-z il^iiy - z)-z i)2(r)(j){t) drdydt. 



(3.54) 



Equation (3.45) now follows from (3.52) and (3.54) since 

lim 0*^(2:)= lim lim (a"''^(z) + '"^(z)). 

From (3.43) and (3.451, and by using Lebesgue's dominated convergence theorem applied to the sequence 
(a'^(z))e>o, we obtain 



lim 



Kg [x^t) ■ ip (x^ — dxdt 

OxR+ ' V £ / 

(y- z,y)u{r,y,t)-z ipi{y~'z)-z '4)2{r)4'{t) drdydtdz 



I Hs(Q) I^P JnxYxM+ 



nxR+ Jhs{o) 



i-pj / a{y - z,y)u{r,y,t)-z ipiiy - z)-zdydz ^p2{r)(l){t)drdt, 



(3.55) 



where we have changed the order of integration in the last step. After shifting the domain of integration 
the inner integral of the right hand side of equation (3.55), we obtain 

o^{y - z,y)u{r,y,t)-z ipi {y- z)-zdy 



ly-z 



^{y,y + z) u{r,y + z,t)-z tpi {y)-zdy 



= / a{y,y + z)u{r,y + z,t)-z ^pi{y)-zdy, 



(3.56) 



where in the last step the fact that the integrand is ^-periodic in y was used. Substituting (3.56 1 in equat 
(3.55), then by changing the order of integration we obtain 



ion 



lim 



axB+ 



Kl,u%x,t)-^/j(x,'^,t^ 



dxdt 



aiy,y + z) -r-r^u{r, y + z,t)dz- ipi {y)dy -02 {r)(j){t)drdt 
nxR+JyJHsio) \z\ 



Bs,iu{r, y, t) ■ 0(r, y, t) drdydt, 



where 



Bs,iu{x,y,t) 



a (y, y) — — ^ ^1^3 — — u{x, y, t)dy 

Hs{y) \y-y\^ 



(3.57) 
(3.58) 



2 

and Kg ^u'^ Bs^iu{x, y, t). 

Next we evaluate the two-scale limit of Kg^^'^- Let ■0 be a test function in J'. Then by using (2.171, 
replacing v(x) with u^(x,i), we obtain 



nxR+ 



KUu' 



dxdt 



7nxR+ JHs(O) 



{x,t) ■ 4) {x, ^,t^ 

fXX \ Z® Z F/ ,\ , f ^ ,\ 

Xn[x + ez)a \^^^ + ^) " {x,t) ■ ip [x, -,tj 



— .t 1 dxdt. 



(3.59) 



Bacim Alali and Robert Lipton 



20 



The right hand side of (3.59), after changing the order of integration, is equal to 

[ b^z)dz. 

where is given by 



\Z\ ./OVB+ \£ £ J V £ / 



■zdxdt. 



(3.60) 



(3.61) 



2.1 



SUPe>o ll"" ||l°=([0,T];Lp(J1)3) 



< oo hence there is a constant 



JnxM+ 

For future reference note that from Proposition 
C such that the sequence b^{z) is bounded above by 

\b'{z)\ < C\z\-\ for e > 0. (3.62) 
As before we approximate xn by a sequence of smooth functions ^„ such that ^„ — > xn in ^foc(^^) ^^'^ write 
^n(^) = TT? / Cn(x + £z)a ( — , — h z\ u'^ (x , t) ■ z ( X , — , t] ■ z dxdt . 

\z\-^ JnxR+ \£ £ J \ £ J 

Next using the fact that (w^)e>o two-scale converges to u{x,y,t), we see that for z € Hs{0), 

\im¥(z)^ lim lim 6f, (z) = -p-pr / a {y,y + z) u{x,y,t)- z ip {x,y,t)- z dxdydt. 

s^O «^ooe^0 \z\-^ JnxYxR+ 



(3.63) 



From (3.591, (3.60) and (3.64 1, and by using Lebesgue's dominated convergence theorem, we obtain 

dxdt 



lim 



1 



{x,t) ■ ip (^x, ^,tj 

o^iy^y + z) u{x, y,t)-z {x,y,t)-z dxdydtdz 



'Hs{a) I^P JnxYxR+ 

By changing the order of integration and then using the change of variables y ^ y + z, we conclude that 

'{x,t) ■ ip {x, ^,tj 



(3.64) 



lim 



OxR+ 



dxdt — 



nxyxR+ 



Bs,2u{x, y,t) -tl! {x, y, t) dxdydt, (3.65) 



where 



Bs,2u{x,y,t) = I a{y,y)— — f] — —dyu{x,y,t), (3.66) 

iHsiy) \y - yr 



and we conclude that -f^l 2"^(^7 ^) ~^ Bs^2u{x,y,t). Equation (3.41) follows on writing Bs = Bs^i — Bs,2- 

The operator p^^Bs is a bounded operator on U^^^lY \ C (JVf) . This follows from bounds on p^^Bs^i 
and p~^Bg 2- Given any w in LP^^.{Y; C{fl)^) an application of Minkowski's inequality to 
\\p^^ Bs,iw{x,y)\\j^p (Y-c{n)3) shows that p^^Bs,i is bounded. The boundedness of p^^i?s,2 easily follows 
from its definition. 

□ 



3.3 Strong Approximation of Local Fields in Heterogeneous Peridynamic Media 

In this section it is shown that a rescaling in the y variable of solution of the two-scale problem delivers a 
strong approximation to the solution u'^{x,t) of the form u{x,y,t). This is stated in the following theorem. 



Theorem 3.12. Letu{x,y,t) be the solution of the two-scale problem given in Theorem 3.10 then 



lim II (a;, i) - u{x, -,i)||LP(n)3 = 0, 



(3.67) 



for every t in [0, T] and | < p < 00. 
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From the perspective of computational mechanics the numerical effort necessary to discretize and solve 
for u{x, y, t) becomes much less expensive than direct numerical simulation for u^{x, t) when the length scale 
of the microstructure s is sufficiently small relative to the computational domain. In view of Theorem |3.12| 
the numerical computation of u{x,y,t) and the subsequent rescaling y = f provides a viable multiscale 
numerical methodology. This topic is pursued in a forthcoming paper. 



Proof. We start by writing the dynamics for the rescaled function u{x, - ,t). Making the substitution y 

1 value prob 



in (3.191 delivers the following initial value problem for u{x, 



£ 



H-,iO) 



u{x + y', t) dy' - u(x, t) d£, 



+P Hf) / xni^ + sz)a{-,- + z)^T^^ (u{x,- + z,t) - u{x,-,t)) dz 



(3.68) 



+ p-i(f)fo(x,f,t), 



with u{x, f ,0) ^ uo{x f ) and dtu{x, f ,0) = vq{x, f ). 

We subtract (3.68) from (2.8) to arrive at the differential equation for the difference e^[x,t) — u'^{x,t) 
u{x,'z,t) given by 



dte'{x,t) = A%e'{x,t) + Ale^x,t) +d''{x,t) 



(3.69) 



with the homogeneous initial conditions e'^(a::, 0) = and dtd^ix^t) = 0. Here the forcing term d^{x,t) is of 
the form d'^{x^ t) = p^^{^) {d% i + c?5 2 + '^l) where 



IE f t \ ^ \ z ® z ( X X \ 

ds 1 = Xn(2; + £z)a — h 2; , ,0 u(x + £z, — h z, t) — uix, — V z,t)] dz 

/ / ^ ^ / X X \ Z Z / , X . , X ,\ 

(1 - Xn[x + £z))a -, - + z . .„ u[x, - + z,t) - u[x, - + z,t)) 

\£ £ / ZP \ £ £ / 



d% 



£ £ 



u{x + ^, 



dz. 



t)- I u{x + £_,y,t)dy I d£, 



(3.70) 
(3.71) 
(3.72) 



The forcing term d^(x, t) is regular and vanishes as e — )• 0, this is stated in the following theorem. 

Theorem 3.13. The forcing term d'^{x,t) belongs to C{[0,T]; LP(rt)^) and the sequence {d'^)£ is uniformly 
bounded for < t < T where 



sup sup \\d%x,t)\\LP(n)3 < 00, for - < p < 00, 
s>Qte[a,T] ^ 



(3.73) 



lim ||c?'^(x, t)\\i^p(Q)a = 0, for all t G [0, T] and - < p < 00. 
E— >o 2 



(3.74) 



We provide the proof of Theore m |3.13| at the end of this section. Since A'^ is a bounded linear operator 
on LP(ri)^ it follows from Theorem 3.13 and Proposition 2.1 that the solution e^(a;,i) is explicitly given by 



e^x^t) 



°° (f _ _^2rl+l 



t'o + 



(3.75) 



Thus 



\\e'{;t)\\L.^nr < / E 

■Jo n 



(t-r) 



2n+l 



(2n+ 1)! 



\{AT\\ ||rf^(-,r)|Up(o)3dr 



< 



' ^smh[VM{t-r)) ||d^(,r)|l 



L=(n)3 dr 



(3.76) 
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where in the second inequahty we have used the fact that is bounded above by a positive constant M > 
independent of e. In view of Theorem |3.13| we can apply the Lebesgue dominated convergence theorem to the 
right most inequahty of (3.76) to conclude that lime_j.o ||e'^(-, i)||^p(j2)3 = and Theorem 3.12 is proved. □ 

We conclude this section by proving Theorem |3.13| The theorem is proved by showing that each compo- 
nent of given by p^'^d^s.i^ Pr^4.2> P^^^^s.s belong to C{[0,T]; LP{n)^) and satisfy ([3J3| and ^J^. We 
begin by showing that pj^dg-^ satisfies (3.731 and (3.74) and that p^^dg^ belongs to C{[0,T]; LP{il)^). In 
what follows we use the basic estimate stated in the following lemma. 

Lemma 3.14. For any subset S of fl and v{x,y,t) in C([0, T]; C(ri)'^)) there exists a fixed integer 

independent of e denoted by L > for which 



\v{x,-,t)\Pdx 



i/p 



< 



sup \v{x', —,t)\P dx 



i/p 



per\ 



Proof. The proof is identical to the arguments used in the estimate (2.24). 



(3.77) 

□ 



We begin by showing that pj^l ^ satisfies ( |3.73[ ) and ( |3.74[ ) and that p^Ml, ^ belongs to C([0, T];LP{n)^) 
Let a — maxy j^'gy p ^ {y)ot{jj ,y') and estimate 



Pe d%^i\\Lv(n) < 



Xn{x + ez) — \u{x + ez, -,t) — u{x, -,t)\ dz dx 
n \JHs{o) \z\ e ^ / J 



< 



< 



HsiO) 



'Hs{0) 



a 

R 

a 



X X 

Xn(x + ez)\u(x + ez, -, t) - u(x, -, t)P dx 
e e 



i/p 



dz 



sup 

n x'en 



{X X ^ P 

Xnix' + sz)\u{x' + ez, -,t) - u{x' , "^Olj 



dx 



i/p 



dz 



Hs{0) 



-Je{z,t) dz, 



where fe{z,t) is given by 

X X 

feiz,t) = \\xn{x' + ez){u{x' +ez, -,t) - u{x' , -,t))\\Li^^(Y-c{n)^Y 



(3.78) 



(3.79) 



Here the second inequality in (3.781 follows from the Minkowski inequality and the last inequality in (3.78) 
follows from Lemma 3.14| Next we show that lim£_j.o \fi,{z,t)\ = 0. To see this write 

9e{y,z,t) = sup {xn{x + ez)\u{x + ez,y,t) -u{x,y,t)\} (3.80) 

and note that 

• ge for almost every y £ Y, t e [0,T], and z e ^^(O), 

• < 9s{y,z,t) <2sup^^^\u{x,y,t)\, 

and \im^^Q\f^{z,t)\ — follows from the Lebesgue dominated convergence theorem since u belongs to 
LP^^{Y;C{nf) for every t e [0,T]. Observe next that 



SUp|/e(z,t)| < '2\\u\\i^P(^Y;C(n)^) - '^\\'^\\c{lQ,T]-Ll„{Y-C(n)3))- 



Hence we apply the Lebegue dominated convergence theorem again to find that 

lim ||p7id| i||ip(o) = 



(3.81) 



(3.82) 
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and application of (3.81) to the last line of (3.781 gives 



sup sup \\p^ d| iIIlp(O) < 

*G[0,T] £>0 



(3.83) 



Given <t < t' <T we apply Minkowski's inequality together with Lemma 3.14 to obtain the estimate 



\\p7 d%iit)-p7 d%,it')\\L.in) < 2a / \z\-Uz \ \\u(t) - uit')\\^,^^^y.^c^^^,y (3.J 



Since u belongs to C^{[0,T]; LP^^{Y;C{n)^)) the estimate ( |3.84[ ) implies that (i|i(i) belongs to 
C{[0,T];LP,,{Y;C{W)). 

Now we discuss the boundedness, continuity and convergence of p ^dg ^. The overall approach to demon- 
strating these properties for p~^dg ^ the same as before. Here we point out that the mechanism that drives 
p~^dg 2 to zero with e is the point wise convergence l — xn{x+ez) — )• for every x e fi. The norm bounds and 
continuity properties of u{x, y, t) are then used as before to establish the continuity properties, boundedness 
and convergence of the sequence {p~^dg 2)e- 

The overall approach to demonstrating properties for the sequence (p^^d^)^ is also the same, however 
there are some distinctions that arise in the proof of convergence. In what follows we outline the proof of 
convergence pointing out that the continuity proof and boun ds ar e established as before. We begin noting 
that u belongs to Qp with | < p < oo hence from Proposition 



3.7 



u{x, ^,t) ^ u{x,y,t), 



and from Proposition 3.5 it follows that for any test function ip{x) S (£7) with ^ + ^ = 1 that 



ip{x)u{x^ —,t)dx 
e 



ip{x) / u{x,y,t)dydx, as £ — > 0. 



We write 



where 



WPe dl\\Lp(^Q) 



\h,{x)\Pdx 



i/p 



K{x) = xn{x + ^)A^ {u{x + e, t) - "(^ + y'' rfy') dt 

We apply ( |3.86p noting that -0(0 = Xo(a; + Ofp belongs to L^' for p' < 3 to find that 



(3.85) 



(3.86) 



(3.87) 



(3.88) 



lim hAx) = 0. 

e-i-O 



(3.89) 



Application of Holder's inequality to the right hand side of (3.88) for p' < 3 gives the upper bound 

\ i/p' / „ . \ i/p 



\he{x)\ < A 



d^ 



H-,(0) 



H-,{0) £ y 



1/p' 



H-,{0) 



ler^ d^ 



Xn{x + 



u{x + £,,y',t)dy' 



i/p 



d^ 



(3.90) 



Applying Lemma 3.77 to the first term on the right hand side of (3.90), Minkowski's inequality to the second 
term followed with Holders inequality delivers the inequality 

|/is(x)|<C||/i||^.^^(y^c7(n)3), (3.91) 



where C is a positive constant independent of e. From (3.89) and (3.91) it now follows from the Lebesgue 
bounded convergence theorem that 

lim|lp7irfi|Up(o) =0. (3.92) 
The continuity and boundedness properties for p^^d\^ follow along lines similar to the previous arguments. 
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4 Homogenized Peridynamics 



The strong approximation u(x, ^,t) admits a natmal decomposition into a continuous macroscopic compo- 
nent and a possibly discontinuous fluctuating component. The macroscopic component u^{x,t) is obtained 
by projecting out the spatial fluctuations and the corrector r{x, - ,t) containing the possibly discontinuous 
fluctuations is given by the remainder, i.e., 



where 



and 



—,t)— [x, t) + r(x, —it), 
e e 



u"{x,t) ^ (u) ^ / u{x,y,t)dy 



r{x, -,t) 
e 



u{x, -,t) 

£ 



i"{x,t). 



The weak convergence expressed by Proposition |3 . 5| gives 

lim -pr-T / u"^ (x , t) dx — \im -r—r [ u(x, — ,t) dx ~ -r— [ u^(x,t)dx, 
1^1 Jv 1^1 7y £ \V\Jv 



and 



lim 



o\V\ 



r{x, —,t)dx = 0. 
V £ 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



It is evident from (4.4) that the macroscopic component tracks the average or upscaled behavior of 
the actual field u^. Conversely the macroscopic or "averaged" observations of the actual deformation 
will track the dynamics of . Thus it is of compelling interest to obtain an explicit evolution equation for 

in order to qualitatively account for observations made at macroscopic length scales. In what follows 
we show that averaging the two-scale peridynamic equations over the y variable delivers a coupled system 
for the macroscopic and microscopic components u^{x,t) and r{x,y,t). This coupling is seen to impart a 
history dependence on the evolution of . We express this memory effect explicitly by eliminating r and 
recovering an integro-differential equation in both space and time for . 

In what follows we set [t) = {-it) and r{t) — r{-, t) and we denote spatial averages of fields v{x, y, t) 
taken over the y variable by {v){t) = Jy v{x, y, t) dy. Let the constant 3x3 matrix K be defined by 



if = A 



d^ 



lh-,[0) W 

and the coupled dynamics for the evolution of {t) and r{t) is given by the following theorem. 
Theorem 4.1. 

{p-^)KLu"{t) + {p-'Bsr){t) - K{p-\){t) + {p-'b){t), 
(p-i - (p-i)) KLu"it) + (p-'Bsrit) - {p-'Bsr)it)) 
K {p-\{t) {p\){t)) + {p-'b{t) {p-'b){t)) , 

(uo), u"iO) = (vo), r(0) = Wo - (uo), and f(0) 



(4.6) 



m 



(4.7) 
(4.8) 



with initial conditions {0} 

Proof. We write u{x, y, t) = {x, t) + r{x, y, t) and substitute this into the two-scale peridynamic equation 
(3.19). Next multiply both sides of ( 3.19[ ) by p^^ and then take the average both sides of (3.19) with respect 
to the y variable. The equation for given by (4.7) follows noting that (r){t) = and 



(f)(t) = 9,2(r)=0. 



(4.9) 



where the operations of differentiation and integration commute since since u S C'^{[0,T]; LP{Y;C{ft)^)). 



The equation (4.8) follows on substitution of (4.7) in (3.19) 



□ 
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Let 



Now we obtain an evolution equation for by eliminating r from the system given by (4.7) and (4.81 
Cr(i) = p-^Bsrit)-{p-'Bsr){t)-K{p-'rit)-{p-'r)it)), 



and (4.8) becomes 



m = Cr{t) + (p-i - (p-i)) KLu"{t) + p-'bit) {p-'b)it). 



(4.10) 
(4.11) 
(4.12) 
(4.13) 

Proceeding as before one finds that C is a linear operator on L'p^j.{Y; C{il)^) and v{t) and w{t) are given 



Since (4.11 1 is linear we set r — v + w where 

v{t)^Cv{t) + {p-^-{p-'))KLu"{t), 



with initial conditions i;(0) = 0, '!j(0) = and 

w{t) = Cw{t) + p'H{t) - {p-^b){t), 
with initial conditions w{Q) ^ uq = uq — (uq) and w(0) — vq — vq ~ {vq 

by 

v{t) = (Vcy^ j\\n\i({t-T)Vc){p~^ -{p-^))KLu"{T)dT 

w(t) — cosh tVC iiQ + (^VC^ sinhtVCwQ 

+ (Vc) ^ f sinh ((i - r)\/c) {p-^b{T) - (/9"^5)(t)) dr. 



Let 



IC = {p-'Bsr){t)-K{p-'b){t), 



(4.14) 



(4.15) 



(4.16) 



then substitution oi r ~ v + w in (4.7) gives the homogenized integro-differential equation for (t) given 
by the following theorem. 

Theorem 4.2. The homogenized deformation (t) is the solution of the integro-differential equation in 
space and time given by 

{p-Y'u"{t) = KLu"{t) + {p-Y^Jc{ycy^ j\mh((t-T)Vc) {p-^ -{p-^))KLu"{T)dT 

+ {p-Y'{f^Ht) + {p-'b){t)), (4.17) 

with the initial conditions {0) — {uq) and 'U^(O) = (wq)- The force generated by the homogenized defor- 
mation f^{t) = f^{-,t) is given by the history dependent constitutive law 

f"{t) - KLu"{t) + {p-^y^JC (Vc) ^* sinh ((t - r)Vc) [p-^ - {p-^)) Klu"{t) dr. (4.18) 

This equation shows that the evolution law for the homogenized deformation is history dependent. 
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